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Abstract 

We study the rj — > 7r°77 decay via the quark-box diagram in the three- 
flavor Nambu-Jona-Lasinio model that includes the £Za(1) breaking ef- 
fect. We find that the 77-meson mass, the rj — > 77 decay width and the 
T] — > 7r°77 decay width are in good agreement with the experimental val- 
ues when the breaking is strong and the flavor 5^(3) singlet-octet 
mixing angle 9 is about zero. The photon energy and the photon invariant 
mass spectra in rj — > 7r°77 are compared with those in the chiral perturba- 
tion theory. 
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Chiral symmetry plays an essential role in the light-flavor QCD. Light pseu- 
doscalar mesons are regarded as the Nambu-Goldstone (NG) bosons of sponta- 
neous symmetry breaking of chiral SU(N)l x SU(N)r symmetry. Properties of 
the NG bosons have been studied in various chiral effective theories successfully. 

Another interesting low-energy symmetry is the axial U(l) symmetry, which 
is explicitly broken by the anomaly. The symmetry breaking is manifested in the 
heavy mass of 7/ meson, which has been studied as the "£7^(1) problem". It is yet 
not clear how strong the anomaly effect is on the pseudoscalar spectrum mainly 
because of complicated interference of the explicit flavor symmetry breaking due 
to the strange quark mass. Two of the authors (M.T. and M.O.) argued in a 
previous paper [|IJ that a strong C/a(1) breaking and consequently a large flavor 
mixing are favorable for the 77 meson. The argument is based on the analysis of 
the rj system as a qq bound state in the three-flavor Nambu-Jona-Lasinio (NJL) 
model. We claim that the 77 mass as well as the 77 — > 77 decay rate supports 
a Ua(1) breaking six-quark interaction much stronger than previously usedP]. 
This interaction causes a large flavor mixing resulting almost pure octet rj ~ rj 8 . 

The purpose of this paper is to extend the previous analysis to another 77 decay 
process, 77 — ► tt°77 and to confirm our claim of the strong {7a(1) breaking. We are 
interested in this process where the internal structure of the rj, it mesons plays 
essential roles because the photon does not couple to the neutral mesons directly. 
Furthermore, chiral perturbation theory (ChPT) gives too small prediction in the 
leading order and higher order terms are expected to be dominant. 

The three-flavor NJL model is one of the phenomenologically successful chiral 
effective models of the low-energy QCD and the lagrangian density we have used 
is as follows, 

C = C + C 4 + C 6 , (1) 
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where ip is the quark field, rh = diag(m u , m^, m s ) is the current quark mass 
matrix, and A a is the flavor U(3) generator (A = ^2/31). The determinant in 
Cq is a 3 x 3 determinant with respect to the flavor indices i,j = u,d,s. The 
model involves the Ul{3) x Ur(3) symmetric four-quark interaction £ 4 and the 
six-quark flavor-determinant interaction £ 6 incorporating effects of the Ua{^) 
anomaly. 

Quark condensates and constituent quark masses are self-consistently deter- 
mined by the gap equations in the mean field approximation. The covariant 
cutoff A is introduced to regularize the divergent integrals. The pseudoscalar 
channel quark-antiquark scattering amplitudes are then calculated in the lad- 
der approximation. From the pole positions of the scattering amplitudes, the 
pseudoscalar meson masses are determined. We define the effective meson-quark 
coupling constants g mq and g nqq by introducing additional vertex lagrangians, 



C vqq = g^ipi^X 1 ^, (5) 

£nqq = 9-Kqq^lb A^^O , (6) 

with A T/ = cos#A 8 — sin#A°. Here is an auxiliary meson field introduced for 
convenience and the effective meson-quark coupling constants are calculated from 
the residues of the gg-scattering amplitudes at the corresponding meson poles. 
Because of the SU(3) symmetry breaking, the flavor A 8 — A components mix with 
each other. Thus we solve the coupled-channel qq scattering problem for the r\ 
meson. The mixing angle 9 is obtained by diagonalization of the gg-scattering 
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Table 1: The parameters and r] — > 7r°77 decay widths for each G 





Gf 


M v [MeV] 


9(M*){deg} 


9mi 


r[eV] 


0.00 


0.73 


138.1 


-54.74 


3.44 


2.88 


0.10 


0.70 


285.3 


-44.61 


3.23 


2.46 


0.20 


0.66 


366.1 


-33.52 


3.12 


2.06 


0.30 


0.63 


419.1 


-23.24 


3.11 


1.71 


0.40 


0.60 


455.0 


-14.98 


3.15 


1.42 


0.50 


0.57 


479.7 


-8.86 


3.20 


1.20 


0.60 


0.54 


497.3 


-4.44 


3.25 


1.04 


0.70 


0.51 


510.0 


-1.25 


3.28 


0.92 


0.80 


0.47 


519.6 


1.09 


3.30 


0.84 


0.90 


0.44 


527.0 


2.84 


3.31 


0.77 


1.00 


0.41 


532.8 


4.17 


3.32 


0.71 


1.10 


0.40 


537.5 


5.21 


3.32 


0.67 


1.20 


0.35 


541.3 


6.02 


3.31 


0.63 


1.30 


0.32 


544.5 


6.66 


3.30 


0.61 


1.40 


0.29 


547.2 


7.17 


3.29 


0.58 


1.50 


0.25 


549.4 


7.57 


3.28 


0.56 


1.60 


0.22 


551.4 


7.90 


3.26 


0.55 



amplitude at the ^-meson pole. The meson decay constant ju (M = tc, K, tj) 
is determined by calculating the quark- ant iquark one-loop graph. The explicit 
expressions are found in [[!]]. 

The parameters of the model are the current quark masses m u = rrid, m s , the 
four-quark coupling constant Gs, the six-quark determinant coupling constant 
Gd and the covariant cutoff A. We take as a free parameter and study r] 
meson properties as functions of Gd- We use the light current quark masses 
m u = rrid = 8.0 MeV (same as in Other parameters, m s , Gd, and A, are 
determined so as to reproduce the isospin averaged observed masses, m^.rriK, 
and f n . 

We obtain m s = 193 MeV, A = 783 MeV, the constituent u, <i-quark mass 
M U) d = 325 MeV and g nqq = 3.44, which are almost independent of Gd- 
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Table 1 summarizes the fitted results of the model parameters and the quan- 
tities necessary for calculating the r\ — > 7r°77 decay width which depend on 
Go- We define dimensionless parameters Crff = —Gd(A/2ti) 4: AN 2 and Gg 1 = 
Gs{A/2ir) 2 N c . When Gff is zero, our lagrangian does not cause the flavor mix- 
ing and therefore the ideal mixing is achieved. The "77" is purely uu + dd and is 
degenerate to the pion in this limit. 

It is found in || that the 77 — ► 77 decay width is reproduced at about C§ = 
0.7. At this value the ratio Gd(ss) /G$ = 0.44 indicates that the contribution 
from Cq to the dynamical mass of the up and down quarks is 44% of that from 
£4. The mixing angle at G c ^ = 0.7 is 9 = —1.3° and that indicates a strong 
OZI violation and a large (u,d)-s mixing. This disagrees with the "standard" 
value 6 ~ —20° obtained in ChPT||. This is due to the stronger {7a (1) breaking 
in the present calculation. The difference mainly comes from the fact that the 
mixing angle in the NJL model depends on q 2 of the qq state and thus reflects 
the internal structure of the rj meson. On the contrary the analyses of ChPT[Q] 
assume an energy-independent mixing angle, i.e., 0(M 2 ) = 9(M 2 ,). 

We are now in a position to study whether the 77 — > 7r°77 decay rate is 
consistent with our picture of 77 with a large OZI mixing. 

The experimental value of the rj — > 7r°77 decay width is 



We evaluate the quark-box diagram given in Fig 1. We follow the evaluation of 
the box-diagram performed in |J. Other possible contributions will be discussed 
later. 

The 77 — > 7T°77 decay amplitude is given by 



{7i°(pMk 1 ,e 1 ) 1 (k 2 ,e 2 )\r ] (p)) = ^vr) 4 ^ + h + k 2 -p)e^ 2 T, u , (8) 




(7) 
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Figure 1: The quark-box diagram for r] — > 7r°77 



where e± and e 2 are the polarization vectors of the photons. T^ v is given by a 
straightforward evaluation of the Feynman diagrams. After calculating traces in 
color and flavor spaces, we obtain 

1 f d A a 6 

= -t^={cose - V2sm6)e 2 g mq g nqq J — — ^U^, (9) 
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(13) 
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Here Ti^ means trace in the Dirac indices and M is the constituent u,d-quark 
mass. Because the loop integration in (9) is not divergent, we do not introduce 
the UV cutoff. Then the gauge invariance is preserved. The inclusion of the 
cutoff that is consistent with the gap equation will break the gauge invariance 
and make the present calculation too complicated. Note that the strange quark 
does not contribute to the loop. 

On the other hand the amplitude T^ v has a general form required by the gauge 



invariance 
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and M n is the rj meson mass. With A and B, the differential decay rate with 
respect to the energies of the two photons is given by 
d 2 T 



dx\dx 2 



M 5 
256tt 2 



+ \\ B \ 2 
4 1 1 



1 2 
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M 2 
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where M n is the ir° meson mass. Though the mass of 77 as a qq bound state 
depends on G c p, we use the experimental value M v = 547 MeV in evaluating 
(18). The Dalitz boundary is given by two conditions: 

1 



M 2 \ M n 
1 1 < X! + x 2 < 1 -. 

t M 2 )- 2 " Mr, 



(19) 



and 



1 / M 2 ' 
Xl +x 2 - 2 Xl x 2 < - \1 - — 



(20) 



In evaluating (10)-(15), one only has to identify the coefficients of p tJ- p u and g^ v . 
Details of the calculation are given in |J . Defining A and B by 



(21) 



we find A and B as 
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We evaluate A and B numerically and further integrate (18) to obtain the r\ — > 
7r°77 decay rate. The results are given in the last column of Table 1 and shown 
in Fig 2. 

Our result is T(r] — > 7r°77) = 0.92 eV at Gff = 0.70 where the rj — > 77 decay 
width is reproduced |1|. At C§ = 1.40 which reproduces the experimental r\ 
meson mass, Y{r\ — > 7r°77) = 0.58 eV. Both are in reasonable agreement with (7). 

In ChPT 0, there is no lowest order 0(p 2 ) contribution to the rj — > 7r°77 process 
because the involved mesons are neutral. Likewise the next order OQo 4 ) tree di- 
agrams do not exist. Thus the 0(p 4 ) one-loop diagrams give the leading term 
in this process, but the contribution is two orders of magnitude smaller than 
the experimental value (7). This is because the pion loop violates the G-parity 
invariance and the kaon loop is also suppressed by the large kaon mass. 

At 0(p 6 ), there exists contribution coming from tree diagrams, one-loops and 
two-loops. The loop contributions are smaller than those from the order 0(p 4 ). 
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Figure 2: Dependence of the r\ — > 7r°77 decay width on the dimension-less cou- 
pling constant G^. The horizontal solid line indicates the experimental value 
T = 0.85 eV and the dashed lines indicate its error widths. 

At 0(p 8 ), more tree diagrams and a new type of loop corrections appear, but the 
loop corrections are also small. 

In @, coupling strengths of the tree diagrams are determined assuming sat- 
uration by meson resonance poles, such as p,u, ao and a-i- This gives 

r(?7 -> vr°77) = 0.42 ± 0.20 eV, (25) 

where the 0.20 eV is the contribution of ao and 02, the sign of which is not known. 
The result is a factor two smaller than the experimental value. 

The contributions of other mesons, such as &i(1235), /ii(1170), foi(1380) || 
and other tree diagrams |J, are found to be small. 

Although these results based on ChPT are not too far from the experimental 
value, it is noted that the higher order 0(p 6 ) terms in the perturbation expan- 
sion are larger than the leading 0(p 4 ) terms and the results contain ambiguous 
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parameters that cannot be determined well from other processes. 

On the other hand in [|l(J the 0(p 6 ) tree diagrams are evaluated by using 
the extended NJL (ENJL) model They calculated three contributions in 
ENJL, namely, the vector and scalar resonance exchange and the quark- loop 
contributions. Their result is T(i] — > 7r°77) ~ 0.5 eV. They further introduced 
the 0(p 8 ) chiral corrections as well as the axialvector and tensor meson exchange 
contributions, and finally obtained T(r] — > 7r°77) = 0.58 ± 0.3 eV. 

The difference between our approach and that in [[K]] follows. The 

ENJL model lagrangian has not only the scalar-pseudoscalar four quark interac- 
tions but also the vector-axialvector four quark interactions. However, the 
breaking is not explicitly included in their model and therefore the 77 — 77' mix- 
ing is introduced by hand with the mixing angle 9 = —20°. We stress that the 
introduction of the £7a(1) breaking interaction is important to understand the 
structure of the 77 meson. 

There is another difference. The coupling constants of the chiral effective me- 
son lagrangian predicted in the ENJL model are parameters of the Green function 
evaluated at zero momenta. On the other hand we evaluate the quantities at the 
pole position of the mesons. 

Calculated spectrum of the photon invariant mass square m 2 ^, for the 77 — > 
7r°77 decay is shown in Fig 3. As this spectrum is compared with those calculated 
by ChPT in ||, we find ours to be similar to the one for d$ = 4.5 x 10~ 2 GeV -2 
in H which involves an additional 0(p 6 ) contribution to the original Lagrangian. 

Spectrum of the photon energy for the 77 — > 7r°77 decay is shown in Fig 4, 
and given in @ in ChPT. Both are also similar, though there is no experimental 
result. 

In our calculation of the 77 — > 7r°77 decay, we evaluate only the quark-box 
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Figure 3: Spectrum of the photon invariant mass m 



77' 



diagram in Fig 1. Since the vector and axialvector four-quark interactions are 
not included in our model, the only other contribution to this process is the 
scalar resonance exchange. In the ENJL model the contribution of the scalar 



resonance exchange is small |10| . We expect that similar result will be obtained 
in our approach. 

If one includes the vector and axialvector four-quark interaction in the NJL 
model, the pseudoscalar meson properties are affected through the pseudoscalar- 
axialvector channel mixing and the model parameters with and without the vector 
and axialvector four-quark interaction are different. We expect that the models 
with and without the vector-axialvector interaction predict similar results for the 
processes involving only the pseudoscalar mesons with energies much below the 
vector meson masses. It is further argued that the contribution of the quark-box 
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Figure 4: Spectrum of the photon energy E^. 

diagram to the 77 — > 7r°7r° process, that is similar to 77 — > 7r°77 , is quite close to 
that of the vector meson exchange in the vector dominance model[|T2||. 

In summary, we have studied the 77 — > 7r°77 decay in the three-flavor NJL 
model that includes the Ua(1) breaking six-quark determinant interaction. The 
77 meson mass, the 77 — > 77 decay width and the rj — ► 7r°77 decay width are 
reproduced well with a rather strong £/a(1) breaking interaction, that makes 
771 — ?7s mixing angle 9 ~ 0°. Since the 77' meson is expected to be sensitive 
to the effects of the Ua(X) anomaly, it is very important to study the 77' meson 
properties. It should be noted, however, that the NJL model does not confine 
quarks. While the NG bosons, tt, K and 77, are strongly bound and therefore can 
be described in the NJL model fairly well, we do not apply our model to the 
heavy mesons such as p,u and 77'. Further study of the breaking and the 
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r)i — 1] 8 mixing will require a calculation including the confinement mechanism. 
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